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ABSTRACT 



Aims. We have observed the deuterated gas in the high-mass star formation region IRAS 05345+3157 at high-angular resolution, in 

order to determine the morphology and the nature of such gas. 

Methods. We have mapped the N2H+ (1-0) line with the Plateau de Bure Interferometer, and the N2D+ (3-2) and N2H+ (3-2) lines 

with the Submillimeter Array. 

Results. We have detected two condensations in NiD^, with masses of ~ 2-3 and ~ 9 Mq and diameters of 0.05 and 0.09 pc, 

respectively. The high deuterium fractionation (~ 0.1) and the line parameters of the N2D^ condensations indicate that they are likely 

low- to intermediate-mass pre-stellar cores, even though other scenarios are possible. 
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1. Introduction 

The initial conditions of the star formation process are still 
matter of debate. Studies have begun to unveil the chemi- 
cal and physical properties of starless low-mass cores on the 
ve rge of forming low-mass sta rs (Kuiper e t al. ,1996 ; Caselli et 
al. l2002al;l20025 Tafalla et aL 120021; |2005) . demonstrating that 
in the dense and cold nuclei of these cores C-bearing molecular 
sp ecies su ch as CO and CS ar e strongly depleted (e.g. Caselli et 
al. l2002bl ; Tafalla et al. |2002|) . while N-bearing molecular ions 
such as N2H+ and N2D+ (Caselli et al. 2002a; 2002b; Crapsi et 
al. 2005) maintain large abundances in the gas phase, and their 
column density ratio reaches values of ~ 0.1 or more, much 
higher than the cosmic [D/H] elemental abundance (~ 10"^, 
Oliveira et al. 2003). 

The characterisation of the earliest stages of the formation 
process of high-mass stars is more difficult than for low-mass 
objects, given their shorter evolutionary timescales, larger dis- 
tances, and strong interaction with their environments. In order 
to check if the chemical properties peculiar of the earliest stages 
of low-mass stars are valid also for high-mass stars, Fontani et 
al. (2006) observed the deuterated gas in 10 sources selected 
from two samples of high-mass protostar candidates (Molinari 
et al. 1996; Sridharan et al. 2002), that are believed to be the 
closest to the earliest stages of the high-mass star formation pro- 
cess. One of these, IRAS 05345+3157 (hereafter 1), stands out 
because of its very interesting characteristics: it is a luminous 
(1.38 X IO-'Lq at a distance of 1.8 kpc, Zhang et al. 2005) young 
stellar object, embedded in a massive (~ 180 Mq) dusty clump 



(Fontani et al. 2006) in which Molinari et al. (2002) revealed a 
complex structure in the molecular gas observed in the HCO^(l- 
0) line. From IRAM-30m spectra of the N2H+ (1-0) and N2D+ 
(2-1) lines, Fontani et al. (2006) have derived an average deu- 
terium fractionation of ~ 0.01, quite close to the values found 
by Crapsi et al. (2005) in low-mass starless cores. These results 
indicate the presence in the source of molecular gas with phys- 
ical conditions similar to those of low-mass starless cores (i.e. 
r ~ 10 K and nu2 ~ lO^cm"^ ). However, the precise location 
and the distribution of this gas, required to understand its na- 
ture, can be determined only through higher angular resolution 
observations. 

In this letter, we present observations of N2H^ and N2D+ to- 
wards 1, obtained with the Submillimeter AiTay (SMA) and the 
Plateau de Bure Interferometer (PdBI), and we report on the de- 
tection of two compact condensations of N2D^ that have chem- 
ical features typical of low-mass pre-stellar core candidates. A 
full report and a more detailed analysis of the data obtained will 
be presented in a forthcoming paper. 

2. Observations and data reduction 

Observations of N2D+ (3-2) (at 231321 MHz) and N2H+ (3-2) 
(at 279511 .7 MH z) towards 1 were carried out with the SMAQ 
(Ho et al. I2OO4I) in the compact configuration on 30 January 
and 21 February 2007, respectively. The correlator was con- 
figured to observe simultaneously the continuum emission and 
several other molecular lines. The phase center was the nom- 
inal position of the sub-mm peak detected with the JCMT 
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(Fontani et al. 2006), namely R.A.(J2000)=05''37"'52.4' and 
Dec.(J2000)=32°00'06", and the local standard of rest veloc- 
ity Vlsr is -18.4 km s"' . For gain calibration, observations of i 
were alternated with the sources 3C111 and J0530+135. 3C279 
and Callisto were used for passband and flux calibration, re- 
spectively. The SMA data were calibrated with the MIR pack- 
age (Qi 2005), and imaged with MIRIAD (Sault et al. Il995l) . 
Channel maps were created with natural weighting, attaining 
a resolution of: 3."7x3."0 for the NaD^ (3-2) channel map; 3. 
"0x2."8 for the 225 GHz continuum image; 2."7x2."0 for the 
N2H+ (3-2) channel map; l."9xl."2 for the 284 GHz continuum 
image. 

We observed the N2H+ (1-0) line at 93173.7725 MHz to- 
wards 1 with the PdBI on 11 and 20 August 2006, in the D con- 
figuration, and on 3 April 2007 in the C configuration. We used 
the same phase reference and Vlsr velocities as for the SMA ob- 
servations. The nearby point sources 0507-1-179 and 0552-1-398 
were used as phase calibrators, while bandpass and flux scale 
were calibrated from observations of 3C345 and MWC349, re- 
spectively. For continuum measurements, we placed two 320 
MHz correlator units in the band when making the observations 
in D configuration, and six 320 MHz correlator units in C config- 
uration. The N2H^ lines were excluded in averaging these units 
to produce the final continuum image (at ~ 96095 MHz). The 
synthesized beam size of the N2H^ channel map was 3."2x3. 
"4, while that of the continuum was 3."lx3."2. We stress that 
the observations of the N2H^ (1-0) and N2D^ (3-2) lines have 
approximately the same angular resolution. The data have been 
reduced with the GILDAS software, developed at IRAM and the 
Observatoire de Grenoble. 

The spectra of both N2H^ and N2D^ obtained from the 
cleaned maps have been analysed with the software CLASS, 
with the method described in Sect. 2.1.1 of Fontani et al. (2006). 



3. Results 

3.1. Structure of the source 

The map of the integrated intensity of the N2D^ (3-2) line, su- 
perimposed on the 225 GHz continuum map, is shown in the left 
panel of Fig. [T] The map indicates that the N2D^ emission arises 
from two molecular condensations: an extended clump located 
~ 10" N-E of the map center and elongated ~ 15" in the N-S di- 
rection, and a compact core, ~ 5" in size, located ~ 15" S-E of 
the map center In the following we will identify these conden- 
sations as N and S, respectively, and they represent the targets 
of the present study. The 225 GHz continuum image shows two 
main compact cores inside the SMA primary beam, one approx- 
imately corresponding to the map center and the other located 
~ 10" S-E of the map center. Another compact source outside 
the interferometer primary beam is located ~ 35" south of the 
map center, and it corresponds to the southern continuum source 
detected in the JCMT image by Fontani et al. (2006, see their 
Fig. A. 3). We will call these cores CI, C2 and C3, respectively. 
None of these cores overlaps with the two N2D^ condensations. 
Some extended emission however is detected at ~ 3cr towards N. 
The right panel of Fig.[T]shows the N2H+ (3-2) integrated emis- 
sion superimposed on the 284 GHz continuum emission, both 
observed with the SMA. At this frequency the continuum source 
CI is resolved into a main peak and a fainter secondary peak 
north of the main one. The N2H+ (3-2) line integrated emission 
is compact and well overlaps with core CI, while it is detected 
only at 3cr towards C2, and undetected towards N and S. The 
detection of N2H+ towards CI, and the non-detection of N2D^ 



at that location, is probably due to the fact that the continuum 
image is dominated by the emission of the warm dust, so that 
in the continuum condensation the temperature is too high for 
significant deuterium fractionation. 

The integrated intensity map of the main group of the hyper- 
fin e (hfs ) components of the N2H^ (1-0) line (see e.g. Caselli et 
al. ll995l) observed with the PdBI is shown in Figure|2] The emis- 
sion is extended and irregular, covering an area of ~ 40" in the 
N-S direction, and ~ 20" in the W-E direction. The main emis- 
sion peak is close to C2, and the shape of the emission indicates 
that N2H^ (1-0) and the 93 GHz continuum are fairly-well over- 
lapping. An extended component in both line and continuum ori- 
ented NE-SW roughly corresponds to condensation N. Another 
compact source, C4, is detected ~ 20" south of the map center 
Such a different distribution of the integrated intensity with re- 
spect to that of the N2H^ (3-2) line is probably due to the fact 
that the extended emission is filtered out differently by the two 
interferometers. To compute the amount of the missing flux, we 
have compared our interferometric spectra with the single-dish 
spectra obtained with the IRAM-30m telescope (see Fontani et 
al. 2006). We have resampled the single-dish spectra of N2H^ 
(1-0) and (3-2) to the same resolution in velocity as the inter- 
ferometric spectra. In the N2H^ (1-0) line, the flux measured by 
the PdBI is ~ 2 times less than that measured with the IRAM- 
30m telescope, while in the N2H^ (3-2) line the flux measured 
with the SMA is only one fifth of that measured with the 30m 
antenna. This indicates that the extended emission is much more 
resolved out in the SMA map with respect to the PdBI map. 

3.2. A/2H+ and N2D^ line parameters and deuterium 
fractionation 

In TablelUwe give the N2H+ (1-0) and N2D+ (3-2) line param- 
eters of the spectra integrated over the 3cr rms level of the N2D^ 
emission in condensations N and S: in Cols. 3 - 7 we list inte- 
grated intensity (j TMBdV), peak velocity (Vlsr), FWHM, opac- 
ity of the main component (r^), and excitation temperature (Tex) 
of the lines. The integrated intensities have been computed over 
the velocity range given in Col. 2, while for the other parameters 
we have adopted the fitting procedure described in Sect. 2. 1 . 1 of 
Fontani et al. (2006). We do not give the line optical depths and 
Tex for the N2D^ lines because in both spectra the uncertainties 
are comparable to the values obtained. Because of this, we have 
fitted the lines forcing the optical depth to be 0. 1 . 

Linewidths and peak velocities of both N2H+ and N2D^ are 
very similar in both condensations, indicating that they are trac- 
ing the same gas. In particular, the lines are ~ 0.8 km s"' broad 
for both molecular species, i.e. nearly half of that observed with 
the IRAM-30m telescope. This can be due either to extended 
components with different velocity that have been resolved out, 
or to a decrease in turbulence going from the pc-scale to the 
sub-pc scale. However, it is interesting to notice that the lines 
are broader than the typical N2H^ and N2D+ lines observed to- 
wards low-mass pre-stellar cores, for which values of ~ 0.2-0.3 
km s ' are found on comparable linear scales (Crapsi et al. 2005, 
Roberts & MiUar 2007). 

We have derived N2H^ and N2D^ total column densities, 
A^(N2H+) and A^(N2D+), following the method described in the 
Appendix of Caselli et al. (: 2002bl) . which assumes a constant 
Tex. A^(N2H+), A^(N2D+) and £)ftac= A^(N2D+)/A?(N2H+), derived 
for N and S, are listed in Cols. 2, 3 and 4 of Table |2] respec- 
tively. These have been obtained using the values of Tex given 
in Table[T] For the N2D^ (3-2) lines, we have assumed the exci- 
tation temperature of the N2H^ (1-0) line because we could not 
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Fig. 1. Left panel: map of the emission of the N2D^ (3-2) Hne integrated between -18.37 and -16.7 km s ', obtained with the SMA 
towards 1 (red contours). The two main condensations are indicated as N and S. Contour levels start from the 3cr rms (~ 0.09 Jy 
beam '), and are in steps of 2 cr. The grey scale represents the 225 GHz continuum (first level is the 3crrms - 0.006 Jy beam '; step 
= 3cr rms). The three compact continuum cores are indicated as CI, C2 and C3. The cross corresponds to the map center (1 nominal 
position) and the dashed circle represents the SMA primary beam at 225 GHz (~ 56"). The ellipse in the bottom left corner shows 
the synthesised beam of the NaD^ image. 

Right panel: map of the N2H^ (3-2) line emission integrated between -19.6 and -16.1 km s ' towards 1 (green contours), observed 
with the SMA (first contour - 3cr rms (~ 0.45 Jy beam"'); step - 2cr rms). The grey scale represents the 284 GHz continuum (first 
level is the 3cr rms = 0.018 Jy beam "'; step - 3cr rms). The dashed circle represents the primary beam at 284 GHz (~ 44"). The 
cross indicates the map center and the ellipse in the bottom left corner is the SMA synthesised beam of the N2H^ image. 

Table 1. N2H^ (1-0) and N2D^ (3-2) line parameters. Between parentheses, the uncertainties of the fitting procedure (see text) are 
given. 



Spectral resolution velocity range jTMBdV Vlsr 

(kms-') (kms-') (Kkms"') (kms"') 



FWHM 
(kms"') 



^ ex 

(K) 



N 
N2H+ (1-0) 
N2D+ (3-2) 

S 
N2H+ (1-0) 
N2D+ (3-2) 



0.25 
0.5 



0.25 
0.5 



-29.9; -12.5 8.7(0.3) 

-18.3; -16.7 0.64(0.03) 



-17.915(0.003) 0.732(0.007) 4.132(0.212) 6.7(0.1) 
-17.91(0.02) 0.88(0.014) O.l(fixed) 6.7' 



-26.5; -9.3 
-18.3; -16.7 



10.3(0.5) 
0.82(0.06) 



-17.283(0.005) 
-17.51(0.04) 



0.73(0.01) 
0.88(0.06) 



6.1(0.4) 6.6(0.1) 

0.1 (fixed) 6.6' 



assumed to compute the N2D^ total column density in Table|2] 



derive a reliable value from the fitting procedure taking into ac- 
count the line hyperfine splitting, as already pointed out. In both 
condensations, Df,ac is 0.1 1, which is comparable to the values of 
Dfiac found in low-mass pre-stellar cores by Crapsi et al. (2005), 
following the same method. These values are also comparable to 
those derived by Pillai et al. (2007) in infrared-dark clouds from 
deuterated ammonia. However, their observations are related to 
the very cold, pc-scale molecular envelope, and not to compact 
sub-pc scale cores. 



3.3. Nature of the A/2D+ condensations 

The main finding of this work is that the N2D^ emission in 
1 is concentrated in two condensations, both of them charac- 
terised by high values of deuterium fractionation. We discuss 
now which is the nature of these condensations. 

We have computed the angular diameters of both condensa- 
tions assuming that the N2D^ integrated intensity profile can be 
fitted with a 2D Gaussian: the geometric mean of the major and 
minor axis resulting from these fits has been corrected for the 



beam size. Then, the linear diameters, L, have been computed 
using a source distance of 1.8 kpc (Zhang et al. 2005). Then, 
the mass of an equivalent homogeneous sphere of diameter L, 
Mn,h+, has been obtained from A^(N2H+) assuming a N2H^ av- 
erage abundance of 1 .5x10"'" (Fontani et al. 2006). L and Mn2h+ 
are listed in Cols. 5 and 6 of Table|2] respectively. For N and S, 
we derive Mn,h+ ~ 8.7 and 2.5M0, and L ~ 0.09 and 0.05 pc, re- 
spectively. These values, together with a Dtvac ^0.1, suggest that 
both cores may be low-mass pre-stellar cores. On the other hand, 
as stated in Sect. |3.2| both N and S show lines broader than those 
typically observed in low-mass pre-stellar cores. We propose 
three possible scenarios about the nature of the two N2D+ con- 
densations: (i) they represent the residual of the extended cold 
cloud in which the cluster of young stellar objects was formed, 
whose physical/chemical conditions have not been altered yet 
by the cluster members; (ii) they are a group of unresolved low- 
mass starless cores; (iii) they are single starless cores, which are 
either low- to intermediate-mass pre-stellar cores, or the seeds 
of future massive forming stars. 
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Fig. 2. Map of the intensity of the N2H^ (1-0) Hne (green contours) integrated between -21.2 and -14.5 km s"', corresponding to 
the main group of the hyperfine components, observed with the PdBI. Levels range from the 3cr rms, which is ~ 0.02 Jy beam ', 
to 0.26 Jy beam ', in steps of 3cr. The grey scale represents the 97 GHz continuum: the levels range from the 3cr rms (4.2x10 "^ Jy 
beam"') to 4x10"'' Jy beam"\ in steps of 3cr. Below the map, the spectra of N2D^ (3-2) and N2H^ (1-0) integrated over the red 
contours, which correspond to the 3cr level of the N2D^ condensations N and S (see Fig.[T]), are shown. The other symbols have the 
same meaning as in left panel of Fig. [1] The continuum source identified in Fig. [T]as C3 is out of the box. 



The first scenario is very unlikely for both N and S, since 
the N2D^ (3-2) line is expected to trace gas denser than ~ 10'' 
cm""*, while the H2 volume densities measured in the pc-scale 
gas associated with high-mass star formation regions typically 
reach values no higher than ~ lO'' - 10^ cm~^. The third 
scenario is the most probable for S, while for N, which shows an 
elongated structure (see Fig. [T) indicative of possible multiple 
low-mass components, the second scenario is also possible. 
In fact, the shape of N resembles that of Oph A, a molecular 
clump with elongated shape, resolved into several low-mass 
starless cores (Andre et al. 2007). If so, the observed broad lines 
could be simply due to unresolved cores with diffent velocity. A 
detailed comparison between N and Oph A will be done in the 
forthcoming paper However, it is also possible that N is, like 
S, a unique condensation, the observed structure being due to 
shaping by a powerful outflow associated with another cluster 
member. In fact, an outflow driven by the massive source CI 
(~ 20 Mq) has been detected observing the CO (2-1) line (this 
line is part of the large dataset obtained with the SMA that 
will be published in the forthcoming paper), and the red lobe is 
indeed detected at the edge of condensation N. If both N and 
S are single cores, what may cause the observed broad lines? 
The larger turbulence could be simply understood as the result 
of a high pressure environment (see e.g. McKee & Tan |2002[) . 
Despite of this, it is not at all clear which fraction of the line 
width is actually broadened by systematic motions, in particular 
infall or accretion of molecular material onto the pre-stellar 
cores. Only high resolution kinematic studies of the large scale 
gas around these condensations will give us clues on the line 
width partition, thus providing important constraints on the 
dynamical evolution of massive star forming regions. 



Table 2. N2H+ and N2D+ total column densities (A^(N2H+) and 
A^(N2D^)), deuterium fractionation (Z)fi-ac), linear diameter (L), 
and mass derived from A^(N2H^) (Mn,h+) for the N2D^ conden- 
sations N and S of Figure [T] The uncertainties computed fol- 
lowing the standard propagation of the errors are given between 
parentheses. 



source A'(N2H+) 
(xl0"cm-2) 



A'(N2D+) 
(xlO'-cm--) 



£»f, 



L 

(pc) 



(Mg) 



N 
S 



1.9(0.6) 
1.5(0.5) 



2.1(0.3) 
1.6(0.3) 



0.11(0.04) 
0.11(0.04) 



0.09 
0.05 



8.7 
2.5 
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